#Get Data

polys_sf <- readRDS('./clean_data/pixel_values/polys_sf.rds') %>% st_make_valid() %>% as_tibble()
val_pred <- readRDS("./clean_data/pixel_values/val_pred.rds") %>% 
  left_join(polys_sf) %>%
  dplyr::select(c(label, .pred_class, poly_num, area)) 
## Joining, by = c("label", "quadrat", "poly_num")

Area by class

#scaleFUN <- function(x) sprintf("%.2f", x)
val_pred %>% group_by(label) %>% summarize(area = sum(area), n = n()) -> area_summary

ggplot(area_summary, label = n) +
  geom_col(mapping = aes(x = label, y = area, fill= label)) +
  # geom_text(nudge_y = 100000, size = 3) +
  scale_fill_manual(values=c("#b7a386", "#5a6233", "#5e4a31"))# +

  #scale_y_continuous(labels = scaleFUN)

Filter Live Veg

val_pred %>% filter(label == 'live vegetation') -> live_vegetation_data
val_pred %>% filter(label == 'dead vegetation') -> dead_vegetation_data
val_pred %>% filter(label == 'sand') -> sand_data

Live Vegetation

Live Vegetation Exponential Cut

Live Vegetation Regular Cut

Sand

Sand Expotential Cut

Sand Regular Cut

Dead Vegetation

Dead Vegetation Expotential Cut

Dead Vegetation Regular Cut